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We present a study of the bend angle distribution of semiflexible polymers of short and interme- 
diate lengths within the wormlike chain model. This enables us to calculate the elastic response 
of a stiff molecule to a bending moment. Our results go beyond the Hookean regime and explore 
the nonlinear elastic behaviour of a single molecule. We present analytical formulae for the bend 
angle distribution and for the moment-angle relation. Our analytical study is compared against nu- 
merical Monte Carlo simulations. The functional forms derived here can be applied to fluorescence 
microscopic studies on actin and DNA. Our results are relevant to recent studies in "kinks" and 
cyciization in short and intermediate length DNA strands. 
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INTRODUCTION 

A classic study in elasticity is the bending of beams 
or rods subject to forces and moments 0. At the cellu- 
lar level there are many macromolecular structures like 
actin and cytoskeletal filaments, which are like beams 
in giving rigidity and structure to the cell. Unlike the 
beams studied by civil engineers, these macromolecular 
beams are subject to thermal fluctuations. The purpose 
of this study is to look at the role of thermal fluctua- 
tions in shaping the elastic properties of macromolecu- 
lar beams. 

We work within the wormlike chain model@,i] , which 
has been known to describe double stranded DNA Q 
as well as actin filaments For clarity we consider an 
experiment in which one end of the molecule is fixed 
at the origin and its tangent vector at the same end 
is constrained to lie along the z direction. We wish 
to know the number of configurations (counted with 
Boltzmann weight) that will result in the final tangent 
vector it. In this paper we present an approximate 
analytical study of this statistical mechanical problem. 

Some earlier treatments of this problem 0, 0] restrict 
to small bending angles, so that the polymer is essen- 
tially straight. One can then replace the sphere of tan- 
gent directions by the tangent plane to this sphere. This 
gives a good account of small bending angles. However, 
there is considerable experimental and theoretical inter- 
est in large bending angles to understand cyciization of 
DNA [8H10|. Ref. developed a new approximation 
technique that works even for large bending angles, say 
of order tt/2. The treatment of is general and in- 
cludes applied forces and torques. In this paper we ap- 



ply the general theory to a special case to illustrate its 
use: we treat the pure bending elasticity of a semiflexi- 
ble polymer not subject to stretching forces or twisting 
torques. 

Two different experimental techniques are possible to 
probe the elastic properties. 

Distribution of tangent vectors: One can tag the ends 
of the molecule with flourescent dye [12] to determine 
the direction of the initial and final tangent vectors. 
By fluorescence video microscopy, one finds the angular 
distribution P(8) of the bending angle 9 , where 6 is 
defined by cos 9 = ti.tf. This experimental technique 
has been used to study actin in two dimensional studies 
in Ref.[l2|- The angular distribution of 9 gives us the 
clastic properties of the molecule and one can compute, 
for instance, the probability for a given bending angle 
9 and compare with the theoretical expectation. 

Measuring moment angle relations: A more inva- 
sive experimental technique is to tether one end of the 
molecule, attach a magnetic bead to the other and ap- 
ply bending moments to the molecule by varying the 
direction of an applied magnetic field. Note that we 
do not constrain the final position of the molecule x(L) 
but only its final tangent vector t/. A uniform magnetic 
field will result in a pure bending moment without ap- 
plying any stretching force. Plotting the bending mo- 
ment vs the bending angle gives another experimental 
probe of the elastic properties. In this paper we derive 
the predictions of the wormlike chain for both these ex- 
perimental situations. 

We first derive an approximate analytical formula for 
the free energy allowing for thermal fluctuations of a 
wormlike chain polymer. We then plot the expected 
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theoretical distribution of bending angles and the theo- 
retically expected moment-angle relation. We conclude 
with a discussion. 



MECHANICS AND FLUCTUATIONS 

In the simplest wormlike chain model, we model the 
polymer by a space curve x(s). x(s) describes the curve, 
and i(s) = ^P, its tangent vector, s is the arc length 
parameter along the curve ranging from to L, the 
contour length of the curve. x(0) = since one end is 
fixed at the origin. The tangent vectors at both ends 
f(0),t(i) are fixed to U and if respectively. 

The mathematical problem we face is to compute the 
partition function 



Q(U,i f ) = J^exp- 
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In Eq.([T]), the sum is over all allowed configurations of 
the polymer, those which satisfy the boundary condi- 
tions for the tangent vector at the two ends: ti = i(0) = 
z,tf = t(L). The energy functional is given by 
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where A is an elastic constant with dimensions of energy 
times length. The quantity L p = A/kT is the persis- 
tence length. For example actin has a persistence length 
of about 16/im and an elastic constant of 6.7 x 10~ 27 in 
units of Nm 2 . Lengths of actin filaments can go up to 
hundreds of microns, so the molecule is of intermediate 
flexibility. 

We will compute the partition function assuming that 
near the stiff limit (L not much larger than L p ), the 
sum over curves is dominated by configurations near 
the minimum of the energy. 

The minimum energy configurations of the polymer 
are those where the tangent vector t(s) describes the 
shortest geodesic connecting ti to tf (Fig. 1). The 
minimum energy is easily computed to be 
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which is a purely mechanical contribution to the 
energy [111]. 

To go beyond mechanics and include the effect of 
thermal fluctuations (Fig. 1), we use an approximation 
described in [ll| . The approximation consists of using 
an expansion of the energy about the minimum energy 
configuration and retaining fluctuation terms about the 




FIG. 1. Two configurations of a molecule with fixed end 
tangent vectors U and tf. (left) Minimum energy configu- 
ration and (right) the same configuration is shown slightly 
perturbed by thermal fluctuations. 



minimum to quadratic order. This calculation was done 
in general form in Eq.(53) of Rcf. 
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for a single 

molecule subject to forces and twisting torques. The 
determinant of the fluctuation operator O is calculated 
in ll|. In the present case, this fluctuation determi- 



nant can be calculated either by specialising 1^.(53) of 
(llj to the case of zero force and torque of [ll| or by 
the method described in [l3| in the context of Brown- 
ian motion on the sphere. As is to be expected, both 
methods give the same answer 
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Performing the Gaussian integration over the 
fluctuations we find the distribution function 
1/^/ (DetO) exp — E r£ n and we arrive at the ap- 
proximate formula 



Q(0) = 
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for the partition function as a function of the final angle 
9, where Af is a normalization constant to be determined 
by the normalization condition 



[ Q(6)sin9d9 = 1 
Jo 
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which includes the measure sin 8 on the sphere. 

The result ([5]) implies that the number of configura- 
tions with a final tangent vector tf making an angle 9 
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defined by arccos (U-tf) is given by P{9) = Q(9)sm9 
leading to the simple analytic approximate formula for 
the bend angle distribution 



P(8) = —V6 sin0exp| 
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for 9. Note that this is a closed analytic form rather 
than an infinite series [l3| . It is thus suitable and con- 
venient for experimental comparison. 

In some earlier works 0, 0] stiff polymers were con- 
sidered as perturbations about the straightline. This 
effectively describes the tangent vector as varying on 
a tangent plane to the unit sphere. This leads to the 
planar formula 
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for the distribution of final angles. Such a planar ap- 
proach neglects the curved geometry of the sphere. One 
also encounters a hybrid formula [H, [U| 
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in which one computes Q{9) using the planar approxi- 
mation and includes the correct curved measure sin 6d6 
on the sphere. The present treatment (unlike the ear- 
lier planar ones) deals with geodesies on the sphere of 
tangents and thus takes into account the curvature of 
the sphere both in the computation of Q{9) and in the 
measure sin (9. In the limit that the bending angle 9 is 
small, Eq.([7]) reduces to the earlier planar and hybrid 
formulae (|8l9p . Our formula for the bend angle distri- 
bution has a wider range of applicability compared to 
the earlier ones: it is valid for contour length compara- 
ble to the persistence length (stiff and the semiflexible 
regime) and also correctly describes the rare events in- 
volving large bend angles. 



keeping L — Nl and L p 
the wormlike chain model. 



fixed, 



one recovers 



A frequency distribution of 
the number of configurations with angle 9 between the 
initial and final tangent vectors was generated. 

Figure 2 shows a comparison between our predicted 
analytical form, the planar formula and the results of 
computer simulation. In the stiff regime (/3 = L/L p less 
than 1) cq. ([7]) is virtually exact, in the semiflexible 
regime (/3 around 5) it remains a reasonable approxi- 
mation. Thus, in an experiment which takes repeated 
snapshots of the angle between the final and initial tan- 
gent vectors of a molecule, we expect to find a distri- 
bution of angles given by Eq. (JT)). This is a prediction 
of the wormlike chain model and it goes beyond earlier 
studies based on classical elasticity [9(. 
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FIG. 2. (Color online) Comparison with simulations of P(9) 
vs 6 for a range of flexibilities. Plotted are the simulation 
data (blue dots), the planar formula (thick red line) and our 
analytic formula Eq. (Q)) (thin black line). Notice that our 
closed form formula gives a reasonable fit to the simulation 
data even for j3 = L/L p as large as 5. 



MOMENT ANGLE RELATIONS 



DISTRIBUTION OF BEND ANGLES 

We have performed Monte-Carlo simulations using 
the Kratky-Porod model for comparison with the an- 
alytical form ([7]). We randomly generated I0 6 config- 
urations of the polymer chain assuming it to be built 
of N identical short straight segments each of length I. 
Given one segment, the direction of the next segment 
was assumed to be uniformly distributed on a cone of 
semivertical angle A9 = .1 radian around the preced- 
ing segment. In the continuum limit when the number 
of segments goes to infinity N — > oo, I — > 0, A9 — > 



The partition function Q(9) can be converted into 
a free energy by the formula (we drop 9 independent 
terms as they are additive constants in the free energy) 

F(9) = -kT\nQ(9) = — - — In — (10) 
2L 2 sin 9 

and used to calculate the bending moment M needed to 
bend the final tangent vector through an angle 9. The 
approximate closed form predicted by our theory is 



Ijr A9 kT 1 

M = ---(-- cot9) 



(11) 
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The first term in (1111) represents the mechanical elastic 
energy and the second (proportional to kT) is due to 
thermal fluctuations. 
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FIG. 3. (Color online) Bending moment vs bending angle 
theoretically expected. The thin blue line shows the me- 
chanical elastic response, and the thick red line includes the 
effect of thermal fluctuations. The moment is plotted in 
units of 10~ 21 N-m. The values of A and kT are chosen to 
reflect actin filaments under physiological conditions. 

As Fig. 3 shows, the mechanical contribution is 
Hookean. Within a purely mechanical approach, when 
the bending angle reaches 7r, the geodesic connecting ti 
to tf is no longer the shortest one and the configuration 
becomes unstable and buckles to a lower energy config- 
uration. This picture is drastically altered by thermal 
fluctuations. As shown by the red curve in Fig. 3, the 
buckling happens at an angle smaller than it. The ther- 
mal fluctuations have the effect of "softening" the linear 
response to external bending. 



clastic rod, such as a poker or a ball point pen refill, 
one sees that if the bending angle exceeds 7T, the con- 
figuration becomes unstable because the tangent vector 
no longer traces the shortest geodesic. The poker then 
buckles into a new configuration with lower energy. As 
shown in Fig. 3, thermal fluctuations cause the buckling 
to set in at bending angles considerably smaller than tt. 

For small bending angles 9, one can Taylor expand 
the free energy (fTOf and find that the effect of thermal 
fluctuations is to effectively reduce the bending elas- 
tic constant A by kTL/6. The fluctuations effectively 
soften the elastic response of the polymer. For larger 
bending angles, one can no longer think of the fluc- 
tuations as simply renormalising the bending modulus 
A, since the form of the moment angle relation is non 
Hookean. 

While we have used actin typical example of 

a semiflexible polymer, our study is also relevant to 
DNA and efforts to understand the bending elasticity 
of small segments. Experiments and all atom simula- 
tions d, @, Q are performed to understand the forma- 
tion of "kinks" and cyclization of short and intermedi- 
ate length DNA strands. Our present approach is an- 
alytical and computes experimentally relevant quanti- 
ties characterising the bending elasticity of semiflexible 
polymers within the wormlike chain model. We expect 
our results to interest researchers studying actin as well 
as DNA and other biopolymers. 

One of us (SS) acknowledges discussions with M. San- 
tosh on all atom simulations of DNA. 



CONCLUSION 

We have analyzed the bending mechanics and fluctu- 
ations of semiflexible polymers. Our central goal in this 
paper has been to bridge the gap between mechanics 
and statistical mechanics by taking into consideration 
thermal fluctuation effects to quadratic order around 
mechanically stable configurations. To summarize, our 
main results are 1) an explicit formula for the bend an- 
gle distribution expected from the wormlike chain model 
and 2) an explicit analytical formula for the moment an- 
gle relation. 

Our calculational technique assumes that there is a 
unique shortest geodesic connecting U and tf. This as- 
sumption breaks down at 9 = tt. Our analytical formula 
([7f should only be used away from the point 9 = tt 
where there is a spurious divergence. For a classical 
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